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We introduce a new dynamic light scattering method, termed photon correlation imaging, which 
enables us to resolve the dynamics of soft matter in space and time. We demonstrate photon corre- 
lation imaging by investigating the slow dynamics of a quasi two-dimensional coarsening foam made 
of highly packed, deformable bubbles and a rigid gel network formed by dilute, attractive colloidal 
particles. We find the dynamics of both systems to be determined by intermittent rearrangement 
events. For the foam, the rearrangements extend over a few bubbles, but a small dynamical corre- 
lation is observed up to macroscopic length scales. For the gel, dynamical correlations extend up to 
the system size. These results indicate that dynamical correlations can be extremely long-ranged in 
jammed systems and point to the key role of mechanical properties in determining their nature. 

PACS numbers: 64.70.pv,82.70.-y,82.70.Rr 



The slow relaxation of molecular fluids, colloidal sus- 
pensions and granular materials becomes increasingly 
heterogeneous 0,110, IIH as the glass or jamming \(\ 
transitions is approached. In these systems, structural re- 
laxation occurs through the correlated motion of clusters 
of particles, leading to spatial correlations of the local 
dynamics. Simulations [3[ and experiments []], 0, 0, H, Q 
of molecular and colloidal supercooled fluids show that 
the range of dynamical correlations is typically moder- 
ate, extending over at most a few particles close to the 
glass transition. Experiments in driven granular sys- 
tems [H|, [Toj| support similar conclusions for athermal sys- 
tems approaching jamming. However, it has been re- 
cently argued [ll| that in jammed materials any local 
rearrangement should propagate elastically over larger 
distances, leading to long-ranged dynamical correlations. 

Testing these ideas is a challenging task. Spatial cor- 
relations of the dynamics can of course be investigated 
in simulations, but the linear size of the simulation box 
is typically limited to 10-20 particles. Optical or confo- 
cal microscopy of colloidal suspensions or imaging of 2- 
dimensional driven granular materials could in principle 
address this issue. In practice, however, the number of 
particles that can be followed is limited. This is because 
the spatial resolution is a fixed fraction of the field of 
view in conventional imaging techniques. Near jamming, 
however, the particle displacements are very restrained, 
such that a high magnification has to be chosen to resolve 
them; this consequently limits the field of view and thus 
the number of probed particles. 

Here, we overcome these limitations by introducing a 
novel optical method that combines features of both dy- 
namic light scattering [l2[ and imaging. The method, 
termed Photon Correlation Imaging (PCIm) [13[, decou- 
ples the length scale over which the dynamics is probed 



from the size of the field of view, thereby allowing long- 
ranged correlations to be measured. We demonstrate 
PCIm by measuring the range, £, of spatial correlations 
in the dynamics of three systems: a diluted suspension 
of Brownian particles, and two jammed materials, a dry 
foam and a colloidal gel. As expected, we find no spatial 
correlations in the dynamics of the Brownian particles. 
For the foam consisting of deformable bubbles that can 
easily slip around each other, £ is of the order of a few 
bubble sizes, but partially correlated motion is also ob- 
served on larger length scales. For the gel consisting of 
rigidly bound hard colloids, the range of correlation is 
comparable to the system size, thousands of times larger 
than the particle size. 

The Brownian particles are polystyrene spheres of ra- 
dius 265 nm, suspended in a mixture of water and glyc- 
erol at a volume fraction cp ~ 10 -4 . The foam is a com- 
mercial shaving foam consisting of tightly packed, poly- 
disperse bubbles (gas fraction « 0.92), a model system 
used to study intermittent dynamics in previous inves- 
tigations [14]. To avoid multiple scattering, we confine 
the foam in a thin cell, of thickness 1-2 bubble layers. 
The foam coarsens with time: the data reported here are 
taken 6.1 h after producing the foam, when the average 
bubble radius is ~ 80 jiva. As an example of a system 
where jamming results from attractive interactions, we 
study a gel made of aggregated polystyrene particles of 
radius 20 nm at (p = 6 x 10 -4 , 32.3 h after the gel was 
formed. Details of the system can be found in [H, [l6| : 
its average dynamics are representative of those of a wide 
class of soft materials. We express all time scales in units 
of r r , the relaxation time of the average dynamics, defined 
by f(T r ) = e -1 , with /(r) the usual average dynamic 
structure factor. For the Brownian particles, foam, and 
gel, r r = 42, 140, and 5000 s, respectively. For all sam- 
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FIG. 1: (Color online) Schematic view of the low-angle, space- 
resolved light scattering apparatus. 



pies, the dynamics are stationary on the time scale of the 
experiments. 

A scheme of the space-resolved dynamic light scatter- 
ing apparatus is shown in Fig.[H The sample cell of thick- 
ness L = 2 or 0.1 mm for colloids and foams, respectively, 
is placed in a temperature-controlled bath (not shown) 
and is illuminated by a collimated laser beam with a di- 
ameter of 8 mm and a wavelength A = 633 nm. A lens 
with a focal length /l = 40 mm is used to form an image 
of the sample onto the detector of a CCD camera, typ- 
ically with magnification 0.8-1.3. We place an annular 
aperture of radius r a = 9 mm in the focal plane of the 
lens, so that only light scattered within a small angular 
interval A6/6 « 0.1 around 6 « 6.4 deg contributes to 
the formation of the image. The images have a speckled 
appearance; each speckle receives light scattered from a 
small volume of lateral extension « A/A# [17[ and depth 
L. Our setup combines features of both imaging and 
scattering methods: similarly to imaging techniques and 
in contrast to usual scattering methods, each location 
on the CCD detector corresponds to a precise location 
in the sample. However, similarly to scattering methods 
and unlike conventional imaging, only light scattered at 
a chosen angle and thus scattering vector q is collected by 
the detector, allowing one to probe the dynamics of the 
sample on a well defined length scale A « q~ x « 1 /im, 
optimized according to the sample features. 

Speckle correlqgraphy [17j and setups similar to the 
one in Fig. [T] [18|, LL9[ have been used in the past to study 
structural properties. Here we focus on dynamical prop- 
erties. We take a time series of CCD images at a fixed 
rate; each image is divided in regions of interest (ROIs), 
typically squares of side 10 — 20 pixels containing ~ 100 
speckles and corresponding to about 50 2 — 150 2 jam 2 in 
the sample. The local dynamics within a given ROI is 
quantified by a two-time degree of correlation [20| : 
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with r the position of the center of the ROI, I p (t) 
the intensity of the p-th pixel of the ROI at time t, 
and (• • -) R oi( r ) tne average taken over all pixels within 
the ROI. The space and time resolved dynamic struc- 
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FIG. 2: (Color online) Dynamic activity maps; the bars cor- 
respond to 2 mm. The color code indicates, in units of stan- 
dard deviations a, the fluctuations of the local cj with respect 
to its temporal average, a): Brownian particles (r = 0.2 s, 
T / Tr = 0.0048); b) foam (r = 3 s, r/r r = 0.021); c) and d): 
colloidal gel (r = 500 s, r/r r = 0.1) during a rearrangement 
event and a quiescent period, respectively. 



ture factor is given by /(t, r; r) = \//? _1 cj(£, r; r), with 
(3 < 1 an instrumental constant [12[ and /(£, r;r) = 
N -1 J2j k exp[iq'(rj(t)-rjfe(t+r))], where the sum is over 
the N scatterers belonging to the sample volume associ- 
ated to the ROI. The usual intensity correlation function 
#2 (t) — 1 measured in traditional light scattering is the av- 
erage of cj(t, r; r) over both t and r. By applying Eq. (pQ) 
to all ROIs, we build a "dynamical activity map" (DAM) 
showing the local degree of correlation, characterizing the 
change in configuration during the time interval [t, t + r\. 

In Fig. [2] we show representative DAMs for the three 
samples; movies displaying sequences of DAMs can be 
found in [21]. While previous works on dynamical het- 
erogeneity have often dealt with the spatial correlation of 
the dynamics on a time scale comparable to the average 
relaxation time r r [5|, |l0|, [22| , here we focus on smaller 
time lags, for which individual rearrangement events are 
best visualized. As a control test, we show in Fig. [2^ a 
DAM for a suspension of Brownian particles. The dy- 
namics exhibit spatial fluctuations, without any obvious 
spatial correlation, consistent with the behavior expected 
for diluted, non-interacting scatterers. A typical DAM 
for the foam exhibiting an intermittent rearrangement 
event is shown in Fig. [2)3. Contrary to the case of the 
Brownian particles, the dynamics is inhomogeneous: the 
bright spot at the top of the field of view corresponds to 
an extended region with a large loss of correlation, while 
the rest of the sample experiences only a marginal change 
in configuration. The DAM movie [2l| reveals that the 
foam dynamics is dominated by similar events. We iden- 
tify these intermittent, localized events with the bubble 
rearrangements driven by the accumulation of internal 
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FIG. 3: (Color online) Instantaneous degree of correlation 
d(£, r, r), normalized by its standard deviation <j(r, r), mea- 
sured at three different locations in the sample, as a function 
of reduced time t/r T (for the sake of clarity, the curves are 
offset along the vertical axis and only a limited portion of the 
data is shown). Data have been corrected for instrumental 
noise as explained in Ref. [23j]. For each sample, curves are 
labeled by the distance in mm with respect to the location 
for the black curve. Top: Brownian particles (the labels 0, 
0.075 and 3.8 refer to the traces from bottom to top) ; mid- 
dle: foam; bottom: colloidal gel. For all samples, r/r r is as 
in Fig. E] 



stress due to coarsening. Spatial maps of the dynamics 
of the colloidal gel are shown in Fig. and d. Simi- 
larly to the foam, the dynamics is due to intermittent 
rearrangement events (c), separated by quiescent periods 
(d). These events, however, extend over strikingly large 
length scales, comparable to the system size, suggesting 
that any local rearrangement must propagate very far 
throughout the network. 

To better evaluate the spatio-temporal correlations of 
the dynamics, we compare in Fig. [3] the time evolution of 
the instantaneous degree of correlation c/(t, r; r) in three 
ROIs, two of which are adjacent, while the third one is 
at least 3.1 mm apart from the first one. Data for the 
Brownian particles are shown in the top panel. The local 
dynamics fluctuates in a noise-like way and fluctuations 
in distinct ROIs appear to be uncorrelated, whatever the 
distance between ROIs. For the foam (middle panel), in- 
dividual rearrangement events such as the one observed 
in Fig. [2)3 result in deep downward spikes, lasting ~ 3 s, 
departing from the average value of cj by several stan- 
dard deviations. The depth of the spikes is a measure of 
the displacement of the scatterers compared to A. The 



cj traces for adjacent ROIs appear to be strongly cor- 
related, while the spikes in the trace of the third ROI 
are almost completely uncorrelated, consistently with the 
typical event size inferred from the DAM of Fig. [2)3. For 
the colloidal gel (bottom panel), we find that essentially 
all ROIs present strongly correlated cj traces, suggesting 
that the ultra-long ranged rearrangement events depicted 
in Fig. [2b are indeed representative of the typical behav- 
ior of the gel. Their duration is of the order of 10 3 s. 

We quantify the spatial correlation of the dynamics by 
introducing a "four point" correlation function G^Ar, r) 
that compares the dynamical activity in two small regions 
separated by Ar. We define G4 as the crosscorrelation of 
the local dynamics 0, H, 
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Here, 5ci(t, r; r) = c/(t, r; r) — (c/(t, r; r)) t and <t(t, r) = 
^ (5ci(t, r; r) 2 )^, with (• • -) t an average over time. Fig- 
ure 2] shows G4 for the Brownian particles (blue squares). 
The spatial correlation decays nearly completely over one 
ROI size, confirming that the dynamics of diluted Brow- 
nian particles are spatially uncorrelated. The spatial cor- 
relation of the dynamics of the foam and the gel are 
shown in Fig. [4] as orange triangles and black circles, 
respectively. Here, a slightly different normalization has 
been chosen: as discussed in Ref. [23|, cj contains a noise 
contribution stemming from the finite number of speck- 
les. The noise of distinct ROIs is uncorrelated and thus 
does not contribute to the numerator of Eq. ([2]); however, 
it does add an extra contribution for Ar = 0, and it con- 
tributes to the a terms in the denominator. Werthus plot 
only data for Ar > and define G4(Ar, r) = bG^Ar, r), 
with b > 1 chosen so that G4(Ar, r) — » 1 for Ar — » 0. 

For the foam, most of the decay of G4 is well de- 
scribed by a slightly compressed exponential, G4 = 
exp[(-Ai?/O p ], with p = 1.27 and f = 0.6 mm (black 
line), corresponding to 8 bubble radii, denoting the size 
of rearrangement events, in agreement with previous es- 
timates [14|. Interestingly, this decay is followed by a 
plateau G4 « 7 x 10 -2 that extends over the full field 
of view, Ar < 4 mm. This indicates that, although the 
spatial correlations of the dynamics are dominated by 
the typical size of rearrangements, the events have also 
a small impact on the rest of the sample. A compari- 
son with younger foams, where the confinement is less 
severe, suggests that this long-range strain propagation 
is due to the strong confinement of the sample investi- 
gated. Elastic strain propagation has a dramatic effect 
in the gel, where the range of dynamical correlation is 
strikingly large: G4 only decays by about 25% over 8 
mm. A compressed exponential fit of G4 (red line) yields 
p = 1.8 and £ = 16.4 mm, demonstrating that the dy- 
namics are strongly correlated over distances far larger 
than any relevant structural length scale: the particle 
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Electron Microscopy (FEM) @ and the PCIm method 
could be adopted to analyze the FEM speckle images to 
monitor the dynamics with spatial resolution. 
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FIG. 4: (Color online) Spatial correlation of the dynamics 
for diluted Brownian particles (G4, blue squares), for a foam 
(G4, orange triangles), and for a colloidal gel (G4, black cir- 
cles). The lines are stretched exponential fits to the decay of 
G4, yielding £ = 0.6 and 16.4 mm for the foam and the gel, 
respectively. For all samples, r/r r is as in Fig. [2] 

diameter is 20 nm and the gel is formed by fractal clus- 
ters of size « 20 /im [15]. The differences observed in 
the spatial correlations of our two jammed systems can 
be understood considering the differences in the response 
to a local change of configuration. The foam consists of 
deformable bubbles, for which the resistance to sliding 
around each other is relatively small. Any local rear- 
rangement is then mostly dissipative, involving the topo- 
logical rearrangement of a few bubbles, and setting only a 
small elastic strain on the surrounding medium. By con- 
trast, in the colloidal gel the rigid connection of particles 
to a stress-bearing backbone entails a long-ranged strain 
field, as both the particles and the bonds between them 
are rigid and do not allow for a localized reorganization. 

Our experiments show that the dynamics of jammed 
systems can exhibit surprisingly long-ranged spatial cor- 
relations, which exceed by far those reported in previ- 
ous works on molecular and colloidal glass formers and 
in driven granular media 0, i, i, i, [|, @, I, E, 0- 
The technique introduced here should provide a flexi- 
ble and valuable tool to characterize spatial correlations 
of the dynamics. PCIm-experiments can be performed 
in any of the classical light scattering geometries, in- 
cluding small angle (where A can be tuned by varying 
ra/fh), wide angle (see Fig. SMI in [Il[), an d diffusing 
wave spectroscopy in the backscattering geometry [24] . 
thereby providing access to a range of length-scales span- 
ning more than 1.5 decades. Ongoing wide angle PCIm 
experiments performed on Laponite [25[, closely packed 
onions [2^, and a fibrin film |2?| suggest that the long 
range dynamical correlations observed here for the gels 
may be very general. In the future, it would be interest- 
ing to compare systematically dynamical correlations in 
jammed (connected) systems and glassy (disconnected) 
systems. Finally, we stress that PCIm could also be ex- 
tended to different kinds of radiation, to probe for in- 
stance molecular glass formers. Indeed, the optical lay- 
out shown Fig. [T] is similar to the one used in Fluctuation 
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